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ABSTRACT 

We calculate the generic spectral signature of an early population of massive stars at high redshifts. 
For metal- free stars with mass above 300 Af©, we find that the combined spectral luminosity per unit 
stellar mass is almost independent of the mass distribution of these stars. To zeroth order, the generic 
spectrum resembles a black-body with an effective temperature of ~ lO^K, making these stars highly 
efficient at ionizing hydrogen and helium. The production rate of ionizing radiation per stellar mass by 
stars more massive than ~ 300Mq is larger by ~ 1 order of magnitude for hydrogen and He I, and by 
~ 2 orders of magnitude for He II, than the emission from a standard initial mass function. This would 
result in unusually strong hydrogen and helium recombination lines from the surrounding interstellar 
medium. It could also alleviate the current difficulty of ionizing the intergalactic medium at z Jj 6 with 
the cosmic star formation rate inferred at somewhat lower redshifts. 

Subject headings: cosmology: theory — early universe — stars: formation — stars: spectra — 
intergalactic medium 



1. INTRODUCTION 

One of the most important challenges in modern cos- 
mology is to understand when and how the cosmic "dark 
ages" ended (Loeb 1998, 1999; Rees 1999). The absence 
of a Gunn-Peterson trough in the spectra of high-redshift 
quasars implies that the universe was reionized at a red- 
shift z > 5.8 (Fan et al. 2000). An early, pregalactic 
generation of stars (the so-called Population HI) has long 
been suspected to control this process (e.g., Carr, Bond, 
& Arnett 1984; Couchman & Rees 1986; Hainian & Loeb 
1997). The feedback of the first generation of stars on 
the ionization history of the intergalactic medium (IGM) 
depends crucially on their initial mass function (IMF). 
Recent numerical simulations presented evidence that the 
primordial IMF might have favored massive stars with a 
mass ^ IO^Mq (Bromm, Coppi, & Larson 1999, 2001; 
Abel, Bryan, & Norman 2000). This prediction relies on 
the existence of a characteristic scale for the density and 
temperature of the primordial gas, and consequently of a 
characteristic Jeans mass of Mj ~ IQ^Mq (see also Lar- 
son 1998). At present, it is not known with any certainty, 
what fraction of a Jeans-unstable clump of gas will even- 
tually be incorporated into the resulting star. This is due 
to the complicated physics of the protostellar feedback on 
a dust-free envelope, and the accretion from it. 

Here we adopt the view that the simulated clump masses 
are indicative of the final stellar masses, and explore the 
implications. Recent observations of the formation process 
of present-day stars seem to support this view. Motte, 



Andre, & Neri (1998) have mapped the nearby p Ophiuchi 
star forming cloud in the 1.3 mm dust continuum. Their 
study identifies a number of dense, gravitationally-bound 
clumps with masses close to stellar values, which appear 
to lack the emission from embedded stellar sources. These 
starless clumps, therefore, might be the direct progenitors 
of individual stars. 

Ultimately, the question of how massive the first stars 
were, has to be answered by observations. It is therefore 
important to ask whether future observatories, such as the 
Next Generation Space Telescope (NGST), could place any 
constraints on the primordial IMF. To address this ques- 
tion, we have constructed models for the interior structure, 
and the emerging spectra of massive Population III stars. 
Since metal-free stars have systematically higher effective 
temperatures [T^ff - 10^ K for M > IOOMq), they are 
very efficient at producing photons capable of ionizing hy- 
drogen and helium. The reionization of the intergalactic 
helium has traditionally been attributed to a population 
of quasars. A generation of very massive Population HI 
stars would provide an important alternative channel, op- 
erating at redshifts higher than what is typically assumed 
for quasars. 

Recently, Tumlinson & ShuU (2000) have discussed the 
case of metal-free stars with a normal (Salpeter-like) IMF, 
up to masses of 9OM0 (see also Ciardi et al. (2000b) for the 
effect of introducing a lower mass cutoff at Mc = 5Mq). 
It is important, however, to extend this mass range into 
the regime of very massive stars with M > IOOM0, since 
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both the nature of the feedback on the IGM, as well as 
the spectral characteristics, change in a significant way for 
this regime. 

In §2, we describe our models for the stellar interior 
and for the spectrum. In §3 and §4 we address the ion- 
ization efficiency and discuss the observational signatures, 
respectively. Finally, §5 summarizes the implications of 
our results. 

2. STELLAR MODELS 
2.1. Stellar Structure 

We construct models for static, metal-free stars in com- 
plete (i.e., dynamical and thermal) equilibrium by solving 
the time- independent equations of stellar structure (e.g., 
Schwarzschild 1958; Kippenhahn & Weigert 1990). The 
integration is performed in the standard way by fitting 
solutions for the core and envelope at a suitably cho- 
sen intermediate point, and the considered mass range is: 
100 < M/Mq < 1000. Stars in this mass range have tra- 
ditionally been termed Very Massive Objects, or VMOs 
(Bond, Arnett, & Carr 1984). For such a massive star 
with the concomitant high interior temperatures, the only 
effective source of opacity is electron scattering, given by 
K ~ 0.2(1 + X) cm^ g~^, where X ~ 0.76 is the mass 
fraction in hydrogen. Other sources of opacity are not 
included in our calculations. 

In the absence of metals, and in particular of the cata- 
lysts necessary for the operation of the CNO cycle, nuclear 
burning proceeds in a non-standard way. At first, hydro- 
gen burning can only occur via the inefficient pp chain. 
To provide the necessary luminosity, the star has to reach 
very high central temperatures (T^ ~ 10^'^ K). These tem- 
peratures are high enough for the simultaneous occurrence 
of helium burning via the triple-a process. After a brief 
initial period of triplc-a burning, a trace amount of heavy 
elements has been formed. Subsequently, the star follows 
the CNO cycle. In constructing our main-sequence mod- 
els, we therefore assume a level of pre-enrichment with 
metals, amounting to a mass fraction oi Z = 10~^. The 
resulting models consist of a convective core, containing 
90-100% of the mass, and a thin radiative envelope. Due 
to the high mass and temperature, the stars are domi- 
nated by radiation pressure, and have luminosities close 
to the Eddington limit: L ~ Xedd — 1-25 x 10'^* erg s~^ 
(M/Mq). 

In this work, we do not include the evolution of the stars 
away from the main sequence, which crucially depends on 
the poorly understood mass-loss mechanism. While this is 
clearly an idealization, it is justified by the fact that the 
stars spend a major fraction of their total lifetime, esti- 
mated to be r ~ 0.007McV-^EDD -- 3 x 10^ yr, near the 
location of the main-sequence. 

In Figure 1, we show the resulting sequence of models 
for Population III stars, and compare to the corresponding 
Population I sequence (calculated with Z ~ 0.02). Our re- 
sults agree well with those of El Fid, Fricke, & Ober (1983) 
in the mass range 100 < M/Mq < 500. It is evident that 
the respective effective temperatures approach an asymp- 
totic value which is independent of mass, T^/f ~ 1.1 x 10^ 
K for Population III, and Tgfj ~ 6.5 x 10^ K for Popula- 
tion I. This characteristic behavior of very massive stars 
derives from the assumptions of hydrostatic equilibrium. 



radiation pressure support with a luminosity close to the 
Eddington limit, and the nuclear physics of hydrogen burn- 
ing. The peculiar behavior of very massive Population III 
stars has already been discussed by a number of previous 
authors (sec Bond, Arnett, & Carr (1984) and references 
therein). Bond, Arnett, & Carr (1984), in particular, have 
pointed out that the effective temperature, Te// ^ 10'"' K, 
is almost independent of mass for M > IOOMq. For com- 
pleteness, we here give a somewhat simplified version of 
the argument. 

The mechanical structure of a massive, radiation- 
pressure dominated star is approximately given by a poly- 
trope of index n = 3 (see Kippenhahn & Weigert 1990). 
If wlz) is the corresponding solution of the Lane-Emden 
equation, and z — {z^/ R)r the dimensionless radial coor- 
dinate, one can express the density and pressure as p = 

PcW^, and P = {ttGR^Pc^ /z^)p'^^^, respectively. Here, 
Z3 — 6.89685 is the dimensionless radius of the polytrope, 
Pc = 54.1825(3M/47ri?^) the central density, and the other 
symbols have their usual meaning. With the polytropic 
value for the central pressure, Pc — ll.lGM^/i?^, one 
can estimate the temperature at the center of the star, 
T, = {3Pc/a)^/\ to be 
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To further constrain Tc, we follow Fowler & Hoyle (1964), 
and consider the global energy balance of the star. The 
average nuclear energy generation rate (in erg s"^ g~^) is 
e — L/M = L edm, and we assume e ~ ecNO- Writing 
e = ec{p/ Pc){T /TcY , with c denoting central values, one 
finds 
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A reasonable fit to the CNO cycle energy generation rate 
close to logT ~ 8.1 is given by 



e~ 1.3 X W^ZXp 
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This form derives from a linear fit to the exact expres- 
sion for ecNO (Kippenhahn & Weigert 1990, p. 165), at 
logT ~ 8.1. Using v — ^ in equation (2) gives e 2± 0.05ec- 
From L ~ ^edd = e^ follows the desired relation 
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Combining equations (1) and (4) gives the mass-radius re- 
lation for very massive stars 



M ~ 370M, 
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from which it can be seen that approximately M oa R^. 
Finally, the effective temperature can be estimated from 

'4 
eff 

the radius, we find 



L 2i -Ledd = 4:TrR^(TT^iy. Using equation (5) to eliminate 
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This analytical estimate for Te//, and the very weak de- 
pendence on mass, are in good agreement with the numer- 
ical results in Figure 1. 

2.2. Emergent Spectrum 

Having determined the basic stellar parameters, we now 
construct model atmospheres and emergent spectra. The 
extreme luminosities very close to the Eddington limit and 
the very high effective temperatures require the sophis- 
ticated effort of using NLTE unified model atmospheres 
with spherical extension and stellar winds, as originally 
introduced by Gabler et al. (f989). This new type of 
model atmospheres provides a smooth transition from the 
subsonic, quasi-hydrostatic photosphere to the supersonic, 
radiation-driven stellar wind, and allows the simultaneous 
treatment of spectral features formed deep in the photo- 
sphere (such as optical continua and line wings) and those 
arising from the outer layers, where the influence of winds 
becomes significant (such as strong absorption edges, far 
infra-red continua and line cores). A detailed discussion is 
given by Kudritzki (1998). 

For our work, we use the new unified model code origi- 
nally developed by Santolaya-Rcy, Puis, & Herrcro (1997), 
and further improved by J. Puis (private communication). 
Bound-bound, bound-free and free-free absorption by hy- 
drogen and helium, as well as Thomson scattering are in- 
cluded as opacity sources. A value of 0.1 has been adopted 
for the number ratio of helium to hydrogen. Changing this 
ratio towards slightly smaller values closer to estimates for 
the primordial helium abundance would not lead to a sig- 
nificant modification of the calculated spectra. 

The atmospheric models are fully specified by the ef- 
fective temperature, the gravity and the stellar radius (all 
taken at the Rosseland optical depth Tross = 2/3) together 
with the mass- loss rate M, the terminal velocity Woo, and 
the parameter (3 which describes the radial slope of the 
velocity field via 

v{r)^v^{l-b/r)'' . (7) 

The constant b is chosen to guarantee a smooth transi- 
tion into the hydrostatic stratification at a prespecified 
outfiow velocity smaller than the isothermal sound speed 
(normally 0.1 Cg). Velocity fields of this form are predicted 
by the theory of radiation driven winds (e.g.. Castor, Ab- 
bott, & Klein 1975; Pauldrach, Puis, & Kudritzki 1986) 
and usually lead to an excellent fit to observed stellar wind 
lines based on radiative transfer and model atmosphere 
calculations (see Kudritzki & Puis 2000). 

At zero metallicity, we expect stellar winds to be of only 
minor importance. However, a radiatively driven wind 
might still arise due to the opacity in the He II resonance 
lines, in combination with the contribution from bound- 
free and free-free absorption. No wind models have been 
constructed so far for this extreme case of zero metallicity. 
Kudritzki (2000) has developed a new algorithm, which al- 
lows to deal with line driven winds at very low metallicity 
down to 10~^ of the solar value. As a rough extrapolation 
from his results we adopt M = 10~^'^Mq yr~^, Woo = 1900 
km s^^ and /3 = 1. This choice of stellar wind parame- 
ters results in an almost hydrostatic, but still spherically 
extended density structure. From a number of additional 
test calculations carried out with different wind parame- 
ters we note that the resulting emergent spectra show only 



a very small parameter dependence as long as the winds 
are weak. For instance, enhancing the mass-loss rate to a 
value of M = 1O~^M0 has little infiuence on the energy 
distribution and line profiles. 

In Figure 2, we present the resulting spectra for differ- 
ent masses. Shown is the specific fiux, Fj^, at the surface 
of the star. At these high effective temperatures, there is 
still a rich spectrum of He II and H lines in the EUV, UV, 
optical, and IR. Some of these lines show emission features 
as a result of an overpopulation of their upper levels due 
to NLTE effects. The He II resonance lines have emission 
wings and central absorption cores, whereas some of the 
He II Pickering lines have a central emission peak on top 
of absorption wings. Such types of line profiles refiect the 
stratification of the ratio of the NLTE departure coeffi- 
cients of the upper to the lower line level and are common 
for very hot stars. A first interpretation of the underlying 
physics has been given in the pioneering paper by Auer & 
Mihalas (1972). Strikingly, the He II absorption edge at 
228 A is very weak or turns into emission for higher masses. 
This is again a distinct NLTE effect which has been dis- 
cussed in detail by Husfeld et al. (1984) and Gabler et al. 
(1989, 1992). 

For M > 300Mq, the spectra become increasingly simi- 
lar. This similarity is even more evident when the resulting 
specific luminosity, L^, =^ AttR^F^, is plotted per unit stel- 
lar mass, as shown in Figure 3. Stars with luminosities 
close to the Eddington limit, L^ ex Xedd oc M, approach 
a universal value for their fiux per unit mass. The pre- 
diction that the total spectral luminosity depends only on 
the total mass of stars but is independent of the initial 
mass function, is a unique feature of a population of very 
massive stars. In §4, we explore the observational conse- 
quences of this result, but discuss first the production of 
ionizing photons from these stars. 

3. FEEDBACK EFFECTS 

The spectrum for a population of very massive, zero- 
metal stars deviates most strongly from the case with a 
Salpeter IMF at short wavelengths, near the peak of the 
corresponding black body at Aem — 250 A. An important 
difference is therefore expected in the rate of producing 
ionizing photons. To quantify this effect for the species 
H I, He I, and He II, we calculate the number of ionizing 
photons per unit stellar mass according to 
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(8) 



where vth is the respective threshold for ionization. We 
here make the somewhat idealized assumption that all 
Population HI stars were more massive than 300Mq. 
Therefore, as demonstrated in the previous section, their 
spectra, when scaled to stellar mass, have a generic form 
which is almost independent of mass. Specifically, we take 
the luminosity per unit mass, L^/M, for a IOOOMq star, as 
given by the dot-dashed line in Fig. 3, to be representative 
for the generic form, and use it in evaluating the integral in 
equation (8) . The resulting production rates are presented 
in Table 1. Comparing these rates to the corresponding 
ones for a Salpeter IMF, as given by Tumlinson & ShuU 
(2000), we find an enhancement in the number of ionizing 



photons by a factor of 10 — 20 for H I and He I, and by 
a factor of ~ 75 for He H. We have verified that using 
LvjM for the SOOMq case (dotted fine in Fig. 3) results 
in ionization efficiencies that are at most a factor of two 
smaUer: 1.2 x 10^^ photons s~^ -^©^ for the ionization of 
H I, 8.6 X 1047 photons s'^ M'^ for He I, and 2.0 x 10''^ 

photons s~^ -^0^ for He II. By comparing to Table 1, one 
can see that our choice of the lOOOM© spectrum as being 
representative for the mass range IOOOM0 > M > SOOMq 
is well justified. 

To illustrate the implications of this result, consider the 
following zeroth-order estimate for the fraction of baryons 
that has to be incorporated into stars, /^, in order to reion- 
ize the respective species. Assuming that ry ionizing pho- 
tons are needed per atom (ry > 1 to account for the effect 
of recombinations) , we find for hydrogen 
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m^TMsNion{^ I) 
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where ryio = (ry/lO), and tms — 2 x 10^ yr is the ap- 
proximate main-sequence lifetime of a very massive star. 
For helium, one has /*(He I) ~ 10^'' 7710 and /*(He II) ~ 
4 X 10~^77io. The corresponding fractions for a Salpeter 
IMF are /. ~ lO-^ryio- 

Consequently, if the first stars are indeed characterized 
by a heavy IMF, the intriguing possibility arises that the 
rare star forming clouds, which originate from the high- 
er peaks in the primordial density field, might have played 
an important role in the reionization of both hydrogen and 
helium at high redshifts. Furthermore, helium might have 
even been reionized before hydrogen. Such a sequence 
of events would be in marked contrast to the standard 
scenarios of reionization (e.g., Gnedin & Ostriker 1997; 
Haiman & Loeb 1997; Gnedin 1998; Ciardi et al. 2000a; 
Haiman, Abel, & Rees 2000; Miralda-Escude, Haehnelt, 
& Rees 2000). Recently, Madau, Haardt, & Rees (1999) 
have pointed out that the known populations of quasars 
and star-forming galaxies are not sufficient to account for 
the required number of ionizing photons at redshifts z ~ 5. 
An early generation of very massive stars could alleviate 
this problem, and contribute a significant part of this un- 
explained deficit. 

4. OBSERVATIONAL SIGNATURE 

4.1. Stellar Spectrum 

Since the first clusters of stars form prior to the epoch 
of reionization, they are embedded within a neutral IGM. 
Consequently, the observed flux shortward of the Lya res- 
onance wavelength, Aq — 1216 A, will be suppressed by 
the neutral IGM before reionization (Gunn & Peterson 
1965) and by the Lya forest after reionization. The sup- 
pression will be complete (due to the neutral IGM) for ob- 
served wavelengths between (1 + Z].oion)-^c and (1 -I- Zs)\a 
and partial (due to the Lya forest) at shorter wavelengths 
(Haiman & Loeb 1999). Here, Ac = 912 A is the Lyman- 
limit wavelength, Zg is the source redshift, and ^reion is 
the reionization redshift. The sharp trough of the ob- 
served flux at {1 + Zs)Xa can be used to unambiguously 
determine the redshift of the source, even in the absence 
of strong spectral lines. To predict the observed spec- 
trum at wavelengths longer than Aq, we consider L{vem), 



the monochromatic luminosity per unit stellar mass at the 
emitted frequency Vem- We again make the assumption 
that the spectrum for a star of IOOOM0 is representative 
for the mass range M > 300Mq. With the flux at the sur- 
face of the IOOOMq star being F{vem), as shown in Figure 
2, one has L{vem) = 47ri?^i^(i^em), where R ~ 13.7i?Q is 
the stellar radius. The observed flux, unattenuated by the 
neutral IGM, at a frequency Vobs — f^em/(l + Zg) is 
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In Figure 4, we show the predicted spectrum for a cluster 
at Zs — 10, assuming a flat universe with Q\ = 1 — Clm = 
0.7 and a Hubble constant of Hq = 65 km s^^ Mpc^^. 
We have taken into account the red damping wing of 
the Gunn-Peterson trough, using the analytical expression 
given by Miralda-Escude (1998). In evaluating the opti- 
cal depth in the vicinity of the Lya resonance, we assume 
•Zroion = 7 and flhh — 0.03, corresponding to a Gunn- 
Peterson optical depth oi to{zs) = 7.4 x 10^. We find that 
the resulting line shape is not very sensitive to the choice 

01 Zreion- 

The peak flux is conveniently located in the spectral 
range between 1 — 5/im, where NGST is expected to reach 
a photometric sensitivity of ~ 0.1 — 1 nJy We plot the 
flux per IO^Mq, but it is straightforward to adjust the 
spectrum to other cases due to the linear scaling with to- 
tal stellar mass. The composite spectrum of a Salpeter 
IMF with stellar masses in the range 0.1 ^ M/Mq <; 100, 
has been worked out by Tumlinson & ShuU (2000), and 
is also shown in Figure 4, normalized to the same total 
stellar mass. In comparing the two spectra, it is evident 
that the observed flux at any given wavelength is a factor 
of 5 — 10 larger for the case calculated with a heavy IMF. 
The difference in absolute flux, however, will not allow to 
discriminate between the two IMFs, since the total stellar 
mass of the observed system is not known a priori. One 
instead has to rely on the difference in colors. 

For this purpose, we define the ratios F12/F23 and 
F23/F45, where F12, F23, and F45 are the integrated fluxes 
between 1 — 2/im, 2 — 3/im, and 4 — 5/im, respectively. 
The two colors can then be expressed as a difference in 
magnitude 

-2.5 log |H , 
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Using the spectra in Figure 4, we find these ratios to be 
F12/F23 — 2.7 and F23/F45 ~ 9.3 in the case of a heavy 
IMF, and F12/F23 ~ 2.5 and F23/F45 c± 7.3 in the case of 
a Salpeter IMF. In terms of magnitudes, one has 



and 
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-1.1 for a heavy IMF 
-1.0 for a Salpeter IMF 



-2.4 for a heavy IMF 
-2.1 for a Salpeter IMF 



Hence, the observed spectrum in the case of a heavy IMF 
is significantly bluer. Redenning by dust is expected to 
be weak for the first star clusters (and is already small 
in some known high-redshift galaxies, see Kudritzki et al. 
2000). Moreover, the two colors differentiate between the 
IMFs in a way that is not degenerate with conventional 
dust. To see this, we ask: Can the presence of dust in the 
cluster with a heavy IMF cause the two colors to simul- 
taneously agree with those predicted in the Salpeter case? 
Applying the average extinction curve for the interstellar 
medium of the Milky Way (e.g., Savage & Mathis 1979), 
we adjust the amount of dust present in the cluster such 
that color A is reddened to agree with the prediction for 
the Salpeter IMF. The corresponding change in color B re- 
sults in AmB — —1.5, compared to Aitib — —2.1 for the 
Salpeter IMF. Unless one invokes rather contrived proper- 
ties for the dust, the photometric signature of a heavy IMF 
therefore cannot masquerade as the one expected from a 
Salpeter IMF. By measuring these colors, NGST might 
thus be able to either confirm or exclude the dominance 
of very massive stars in the early universe. 

4.2. Strong Emission of Recombination Lines 

The hard UV emission from a star cluster with a heavy 
IMF is likely to be reprocessed by the surrounding in- 
terstellar medium and produce very strong recombination 
lines of hydrogen and helium. With Nion as derived in 
§3, the emitted luminosity if™g per unit stellar mass in a 
particular recombination line can be estimated as follows 



cluster at the line wavelengths is given by 
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where pf™g is the probability that a recombination would 
lead to the emission of a photon in the corresponding line, 
u is the frequency of the line and p%lnt and pf?^g are the 
escape probabilities for the ionizing photons and the line 
photons, respectively. It is natural to assume that the 
stellar cluster is surrounded by an ionization bounded H 
II region, and hence Pcont is close to zero (Wood & Loeb 
2000; Ricotti & Shull 2000). In addition, pf^^^ is likely to 
be close to unity in the H II region, due to the lack of dust 
in the ambient metal-free gas. Although the emitted line 
photons may be scattered by neutral gas, they will diffuse 
out to the observer and eventually survive if the gas is 
dust free. Thus, for simplicity, we adopt a value of unity 

To be concrete, let us consider case B recombination 
which yields pf™g of about 0.65 and 0.47 for the Lya and 
He II 1640 A lines, respectively. In arriving at these num- 
bers, we assume an electron temperature of ~ 3 x lO^K, 
and an electron density of ^ 10^ — 10^ cm~^ inside the H II 
region (see Storey & Hummer 1995). With Nion from Ta- 
ble 1, we then obtain Lf™^ = 1.7 x 10^^ and 2.2 x lO^^ erg 
s~^Mq for the recombination luminosity of Lya and He 
II 1640 A per embedded stellar mass. A cluster of 1O^M0 
in stars would then produce 4.4 and 0.6 x\{fLQ in the 
Lya and He II 1640 A lines. Similarly high luminosities 
would be produced in other recombination lines at longer 
wavelengths, such as He II 4686 A and Ha. 

The rest-frame equivalent width of these emission lines 
measured against the stellar continuum of the embedded 
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where Lx is the spectral luminosity per unit wavelength of 
the stars at the line resonance. From Figure 3, we obtain a 
spectral luminosity per unit frequency L^, = 2.7 x 10^^ and 
1.8 x 10^^ erg s^^ Hz^^Mq^ at the corresponding wave- 
lengths. Converting to Lx, we obtain rest-frame equiva- 
lent widths of Wx = SlOOAand llOOAfor Lya and He 
II 1640 A, respectively. These extreme emission equiva- 
lent widths are more than an order of magnitude larger 
than the expectation for a normal cluster of hot stars with 
the same total mass and a Salpeter IMF under the same 
assumptions concerning the escape probabilities and re- 
combination (e.g., Kudritzki et al. 2000). The equivalent 
widths are, of course, larger by a factor of (l+Zg) in the ob- 
server frame. Extremely strong recombination lines, such 
as Lya and He II 1640 A, are therefore expected to be an 
additional spectral signature that is unique to very mas- 
sive stars in the early universe. Tumlinson & Shull (2000) 
have discussed the possible observational role of the He II 
1640 A line in the context of a standard Salpeter IMF. 

The He II 1640 A line would reach the observer un- 
affected by the intervening IGM, since its wavelength is 
longer than that of the Lya transition which dominates the 
IGM opacity. However, the Lya line emitted by the source 
would be resonantly scattered by the intergalactic H I. At 
redshifts lower than the reionization redshift, most of the 
IGM is ionized and only the blue wing of the line would be 
suppressed by the Lya forest. But if the star cluster lies 
beyond the reionization redshift, the intervening neutral 
IGM acts as fog and obscures the view of the source at the 
Lya wavelength. The line photons emitted by the source 
scatter until they eventually redshift out of resonance and 
escape due to the Hubble expansion of the surrounding 
intergalactic H I. As a result, the source creates a faint 
Lya halo on the sky (Loeb & Rybicki 1999). For a source 
at Zs ~ 10, the halo occupies an angular radius of ~ 15" 
and corresponds to a line that is broadened and redshifted 
by ^ 10'^ km s~^ relative to the source. The scattered 
photons are also highly polarized (Rybicki & Loeb 1999). 
A heavy IMF would make the detection of diffuse Lya 
halos around the first star clusters more accessible to fu- 
ture telescopes, and would provide an important tool for 
probing the neutral IGM prior to the epoch of reioniza- 
tion. Observations which focus directly on the bright cen- 
tral source are more straightforward, and should show a 
strong suppression of the Lya line by the damping-wing 
of the "Gunn-Peterson trough" , caused by the IGM scat- 
tering (Miralda-Escude 1998). The very brightest sources 
might be an exception, in that the H II region they gener- 
ate in the surrounding IGM might shift the Lya damping 
wing and allow a portion of the red wing of the emission 
line to be observed (Cen & Haiman 2000; Madau & Rees 
2000). 

5. SUMMARY AND CONCLUSIONS 

We have explored the consequences of a possible early 
generation of very massive, zero-metal stars. These stars 
are almost fully convective, have luminosities close to the 
Eddington hmit, and effective temperatures, Te// ~ lO^K, 



which are almost independent of steUar mass. Using a 
spherically extended, NLTE stellar atmosphere code, we 
have subsequently constructed the emerging spectra. If 
normalized by the mass of the star, the spectrum has a 
generic form, almost independent of mass for M > SOOM©. 
Based on the photometric colors of the composite spec- 
trum from a cluster of Population III stars, we find that it 
would be feasible observationally to constrain the primor- 
dial IMF with future infrared telescopes such as NGST 
(Fig. 4). ^ ^ 

We also investigated the ionizing photon production 
from these massive stars. Compared to a population of 
stars with a normal, Salpetcr-like IMF, the rate of pro- 
ducing ionizing photons is enhanced by an order of mag- 
nitude for H and He I, and by a factor of ^ 100 for He 
II. As a result, the interstellar medium of a star cluster 
with a heavy IMF is likely to generate extremely strong 
recombination lines, such as Lya and He II 1640 A, with 
equivalent widths that are larger by at least an order of 
magnitude compared to a standard IMF. 

Very massive Population III stars might have played a 
significant role in the reionization of the hydrogen and he- 
lium in the IGM at high redshifts. The ability to efficiently 
ionize He II provides an alternative channel for a process 
which has otherwise to rely on the presence of quasars with 
their non-thermal spectra. Equally important is the pos- 
sible contribution of massive stars to the ionizing photon 
budget at z '--^ 5, where the known quasars and galaxies 



are unable to maintain the IGM in an ionized state as long 
as the IMF is assumed to have the standard form (Madau 
et al. 1999). 

The massive stars considered in this work are likely to 
end up as black holes (Fryer, Woosley, & Heger 2000). If 
7-ray bursts (GRBs) originate from the collapse of these 
massive stars, then a heavy IMF will produce a larger 
abundance of GRBs at high redshifts than previously es- 
timated (Ciardi & Loeb 2000). The bright afterglows of 
these GRBs would provide an important probe of the IGM 
at high- redshifts (Lamb & Reichart 2000). 

A valuable contribution to the early stages of this work 
was made by Pierre Demarque, who adapted the Yale 
stellar structure code to compute models of very mas- 
sive zero-metallicity stars; although the results reported 
in this paper were obtained with a more efficient code 
kindly supplied by Dimitar Sasselov, the models of Pierre 
Demarque provided an important check on their correct- 
ness. The authors wish to thank Joachim Puis for the as- 
sistance in using his unified model atmosphere code, and 
Michael ShuU and Jason Tumlinson for making available 
to us their composite Population III spectrum. We are 
grateful to Martin Haehnelt and Martin Rees for help- 
ful discussions. VB thanks Paolo Coppi for financial sup- 
port under NASA grant NAG5-7074, and AL acknowl- 
edges support by NASA grant NAG 5-7768 and NSF grant 
AST-0071019. 
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Table 1 
Production of Ionizing Photons " 



E N,on Y f^/mo 

■0 



[cV] [photons s^i M ^ 



H I 13.6 1.6 X 10"** 

He I 24.6 1.1 X lO"* 

He II 54.4 3.8 X 10""" 



16 


1 X 10""' 


14 


1 X 10-"' 


75 


4 X 10-5 



''The production rates arc calculated for a IOOOMq star. Within a factor of two, these rates apply 
to all stars with mass above 3OOM0. 

Note. — E is the ionization potential, Nion is the production rate of ionizing photons per unit 
stellar mass, Y is the ratio of Nion for the cases of a heavy IMF and a Salpeter IMF, and /* is the 
required baryon fraction that needs to be converted into stars in order to produce 10r;io ionizing 
photons per species particle (H I, He I, or He II). 
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Fig. 1. — Zero-age main sequence of very massive stars. Left solid line: Population III zero-age main sequence (ZAMS). Right solid line: 
Population I ZAMS. In each case, stellar luminosity (in Lq) is plotted vs. effective temperature (in K). Diamond- shaped symbols: Stellar 
masses along the sequence, from IOOMq (bottom) to IOOOMq (top) in increments of IOOMq. The Population III ZAMS is systematically 
shifted to higher effective temperature, with a value of ^ 10^ K which is approximately independent of mass. The luminosities, on the other 
hand, are almost identical in the two cases. 
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Fig. 2. — Spectra of very massive Population III stars. Solid lines: Flux at the surface of the star (in units of erg s ^ cm •^ Hz ^) vs. 
wavelength (in A). Dashed lines: Corresponding black-body spectrum, evaluated at the respective effective temperatures. Lines are due to 
transitions of H and He II in the EUV, UV, optical, and IR bands. The spectra for M = 100, 300, 500, and IOOOMq are remarkably similar, 
differing only slightly in the strengths of the He II and H I ionization edges at 228 A and 912 A, respectively. 
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Fig. 3. — The normalized spectral energy distribution in the continuum. Luminosity per unit stellar mass (in erg s ^ Hz ^ solar mass -') 
vs. wavelength (in A). Solid line: WOMq. Dotted line: 300Mq. Dashed line: 500Mq. Dot-dashed line: IOOOMq. When plotted per unit 
mass, the spectra admit an almost universal form for M > SOOMq. 
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Fig. 4. — The predicted flux from a Population III star cluster at Zs = 10. Observed flux (in njy per 10® Mq of stars) vs. wavelength (in 
/im). We assume a flat universe with H^ = 0.7, and a Hubble constant of Hq = 65 km s~^ Mpc~^. Solid line: The case of a heavy IMF. 
Dotted line: The comparison case of a standard Salpeter IMF, where the composite spectrum is taken from Tumlinson & ShuU (2000). The 
cutoff below \obs = 1216 A (1 + Zs) = 1.34^m is due to complete Gunn-Pcterson absorption. It can be seen that for the same amount of total 
stellar mass, the observable flux is larger by an order of magnitude for stars which are distributed according to a heavy IMF. 



